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Abstract 

The reformulation of the mode-couphng theory (MCT) of the hquid-glass 
transition which incorporates the element of metastability is applied to the 
hard-sphere system. It is shown that the glass transition in this system is not 
a sharp one at the special value of the density or the packing fraction, which 
is in contrast to the prediction by the conventional MCT. Instead we find 
that the slowing down of the dynamics occurs over a range of values of the 
packing fraction. Consequently, the exponents governing the sequence of time 
relaxations in the intermediate time regime are given as functions of packing 
fraction with one additional parameter which describes the overall scale of 
the metastable potential energy for defects in the hard-sphere system. Impli- 
cations of the present model on the recent experiments on colloidal systems 

are also discussed. 
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I. INTRODUCTION 



In Ref. |jl|J^, the notion of metastability was used in the reformulation of the mode 
couphng theory (MCT) of the hquid-glass transition. It was shown that the temperature 
dependence of the transition is smooth, as observed in recent experiments ^j^, showing no 
evidence for the special transition temperature conventionally assumed in MCT. In this 
paper, we apply the basic picture obtained in Ref. to a more reahstic situation by 
considering the effects of the structure of supercooled liquids on the slowing down of the 
dynamics. For a simple system like hard-sphere fluids considered here, we find that the 
observed slowing down of the dynamics occurs over a range of values of the packing fraction. 
This is in contrast to the prediction by the conventional MCT that there exists a critical 
packing fraction for the transition. Consequently, we find that the exponent parameters 
describing the relaxation sequence in the intermediate time regime are given as smooth 
functions of packing fraction with one additional parameter, which is related to the existence 
of metastable defects in the system. 

Much of the recent theoretical progress on the liquid-glass transition problem has been 
inspired by the MCT 0,^. It was particularly successful in explaining the very elaborate 
sequence of time relaxations which has been observed in many experiments [§,0. The con- 
ventional interpretation of this theory assumes a sharp transition as the system approaches 
the critical temperature or the critical density. Although the sharp nature of the transition 



is smeared by the cutoff effect [0, the notion of a critical temperature or density at which 
the dynamics is arrested still remains in the theory ||lT|. For example, the exponents gov- 
erning the time relaxations are defined only at that particular temperature or density. This 
is, however, incompatible with the recent experimental results showing a more universal 
behavior and smooth temperature dependences. Other recent experiments also indicate 
that the existence of the well-defined transition temperature is very doubtful. 

The basic ingredient in the reformulation of MCT in Ref. [|l],0 is the introduction of the 
metastable defect variable and its coupling to the density fiuctuations in the glass transi- 
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tion problem. The defect variable is supposed to be governed by a very long time scale 
compared to the mass density and the metastability is incorporated through a double-well 
type potential energy for the defect variable. By investigating the conditions that the defect 
variable must satisfy in order to be on the critical surface associated with MCT [|12|, the 
authors of Ref. were able to obtain the observed slowing down without having to adjust 
the temperature. The condition for slowing down is given instead by a self-consistent con- 
dition between the low activation barrier for defects and its weak coupling to the density 
fluctuation [HJ^. The temperature dependence of the parameters describing the metastable 
potential for defects turns out to be smooth. The exponents for the time relaxation sequences 
can be determined for a range of temperatures showing a smooth temperature dependence 
compatible with recent experiments. 

In the analysis of those models, however, the wavenumber dependences were suppressed 
as a simplifying assumption. The static structure factor for the density fluctuation S{q) 
was treated in a very simple manner by assuming that it is constant for g < A and zero 
otherwise, where A is the large wavenumber cutoff. In this paper, we generalize the previous 
model by considering a model where the structure factor has a more realistic wavenumber 
dependence. We focus here on the dynamics of hard-sphere fluids for which an accurate 
approximate analytical treatment of S{q) is available. We restrict our analysis here to the 
intermediate time regime before the primary relaxation where much of the analysis can be 
carried out analytically. A detailed MCT analysis has been carried out in Ref. [1^ for this 



time regime, which is very useful in our analysis here. We note that in the previous model 
imj^ the temperature was represented by the parameter ^ which is proportional to A^ as well 
as to the temperature. The appearance of the large wavenumber cutoff in the temperature 
parameter ^ is directly related to the suppression of the wavenumber dependence in the 
model. When one considers the wavenumber dependence of S{q), however, we find that the 
role played by ^ in the previous case is replaced, in the case of hard-sphere fluids, by the 
more physical parameter, the packing fraction. 

The model in Ref. |jl|J^ is based on the fluctuating nonlinear hydrodynamics (FNH) for 



a set of slow variables consisting of the mass density p(x), the momentum density g(x) 
and the defect density nrf(x). As discussed in Ref. [0, the wavenumber dependence can 
be incorporated into the FNH formulation through the terms in the effective Hamiltonian 
F that contains spatial derivatives of the variables. More generally, one can include terms 
where the coefficients depend explicitly on the position x. We construct F such that the 
terms that depend on the density fluctuation 5p yield the correct static structure factor S{q) 
for the hard-sphere system, where 6p = p — po with the average mass density po- In the 
presence of the defect variable, we must also consider the effect of the defect structure on the 
dynamical slowing down. Unfortunately, a microscopic description for the structure factors 
in the presence of defects is not available even for a simple system like a hard-sphere fluid. 
We assume here that the local behavior of defects is similar to that of particles. Thus we 
assume that any possible effect of the defect structure is reflected in the coupling term in F 



between the density fluctuation and the defect fluctuation Sud, where 5nd ^ nd — Ud [0 with 
the metastable defect density rfd- This assumption is reasonable since a vacancy is created 
whenever a particle moves to a different position to make an interstitial. In this paper we 
assume, as a first approximation, that the coupling term in F between 6p and 6nd has the 
same spatial form as the term quadratic in 6p. As we shall see below, this corresponds to 
choosing the static defect-density correlation function to be proportional to s(g) — 1, where 
s{q) = S{q) / {mpo) is the dimensionless form of the structure factor S{q) with m the mass 
of the particle. This quantity is just the static structure factor for the density fiuctuation 
excluding the self-correlation effect. The potential energy for is chosen to have the same 
double well potential as in Ref. [0,0 . The basic picture of Ref. [||,|| is then retained in the 
sense that the slowing down of the dynamics corresponds to the limit of the low activation 
barrier for defects and the weak coupling. 

As mentioned above, MCT successfully explains the observed sequence of time relax- 
ations. This sequence can be summarized as follows: After the microscopic time scale, the 
system goes into the intermediate time regime where the density auto-correlation function 
exhibits the power-law decay, ~ / + Ait~'^ followed by the von-Schweidler relaxation. 



(f) ^ f — A2t^ . The system then crosses over to the primary (a-) relaxation regime which 
can be described by a stretched exponential decay (p ~ exp{(t/rQ,)^}. 



In the intermediate time regime, a detailed analysis is available [|I^ for the general 
MCT equations including the wavenumber dependence. It is shown there that all micro- 
scopic structural details are summarized into one exponent parameter A, which gives the 
power-law and the von-Schweidler exponents a and h. Thus, given a microscopic model, the 
corresponding A can be calculated. In the present case, we find, assuming the system is 
on the critical surface, that A depends upon the parameter describing the overall scale of 
the metastable potential for defects as well as on the packing fraction. This result that the 
exponent parameter A is given as a function of the packing fraction is conceptually very dif- 
ferent from the results of the earlier MCT calculations on the hard-sphere system |17|J16|-[18 



where only the simplest contribution from the density fiuctuation to the dynamic viscosity 
is taken into account. In those cases, the transition is controlled by the packing fraction 
and thus A is defined only at the critical packing fraction just as in a second order phase 
transition. According to the picture presented here, the transition can take place over a 
range of packing fractions with different values of the exponents. 

In Sec. fT|, we construct the FNH model including the wavenumber dependence and 
present the relevant results to our discussion of the slowing down of the dynamics. Since 
the calculation follows closely those in Ref. [jl],|2| for the wavenumber independent model, we 



emphasize here only the important results. In Sec. |T|, we present the results of numerical 
analysis of the model. In Sec. we conclude with a brief discussion. 



II. MODEL 

The FNH approach developed in Ref. is based on the generalized Langevin equations 
for the slow variables, {p, g, n^}. As mentioned in the previous section, the wavenumber 
dependence can be implemented through general spatially-varying coefficients of terms in 
the effective Hamiltonian F[p, g, nj, which governs the static property. In the present model, 
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F is given by 

g, na] = Fk[Sp, g] + F^[6p] + F,[6p, n,] + F,H. (2.1) 
Here F^ is the usual kinetic energy term Fft-[p, g] = / rf^x g^/(2p). For the terms that are 



dependent on the density fluctuation only, we take the Ramakrishnan-Yussouff form [IS 



F4Sp] = ^ jd"^ [p log ^-5p]-^ j rf^xrf^x' 5p(x)C(x - x')5p(x'), (2.2) 
Tfi J po Zvn J 



where ksT is the temperature, m is the mass and C{x) is the Ornstein-Zernike direct 
correlation function for the density fluctuation [^. For the coupling term, we have in 
general 

Fc[5p, rid] = I d^^d'^' n,(x)5(x - x')5p(x') (2.3) 

for some function -B(x) to be specified later. We write the Fourier transform of -B(x) as 
B{q) = J (i^xexp(iq ■ x)i?(x) = B ■ b{q), where the constant B represents a scale of the 
coupling energy with the appropriate dimension for B{q) and the dimensionless 6(g) contains 
the wavenumber dependence of B{q). In Eq. (|2 . 1|) , the potential energy part Fylna] for 
is given by the same double well potential h{nd) with the metastable defect density na as 
in Ref. F^lrid] = J d^yi h^nd), where h'{nd) = 0. We define for later use p = n^h"{nd) 

and u = n\h"'{nfi) which are parameters describing the shape of h{nd). We note that, if 
p 0, the potential h{nd) develops an inflection point at n = n^, which corresponds to the 
potential with very shallow metastable wells and low activation barriers for the defects. 
The static correlation functions can readily be found from Eq. ( ^.1| ): 

(5p(q)5p(-q)) ^ S{q) = ^^^(^ (2.4a) 
(5p(q)5n(-q)) = -^^S{q)b{q) (2.4b) 
(5n(q)5n(-q)) = (2.4c) 

where hq = pq/tti is the average number density. The first equation reproduces the correct 
relation between the static structure factor s{q) and the direct correlation function C{q). In 
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calculating Eq. ( p.4| ), we neglected the terms of order B'^ / jj, compared to those of order 
which is consistent with the self-consistent condition obtained in Ref. |I| for the slowing 
down of the dynamics: The slowing down corresponds to a weak coupling [B —>■ 0) and a 
shallow metastable well for defects (/i —>■ 0) with the condition B As indicated in 

Sec. |, the coupling energy term is constructed such that the static defect-density correlation 
function behaves like s{q) — 1. This is achieved in Eq. (p.4| ) by taking 

h{q) = -noC{q). (2.5) 



As one can see from Eqs. ( p.2|) , (|2.3|) and ( p.5|) , the coupling term Fc has the same form as 
the quadratic interaction term between the density fluctuations. 

The construction of the generalized Langevin equations for a general set of slow variables 



is described in Ref. |^T|. The main ingredients in the construction are the effective Hamilto- 
nian and the Poisson brackets between the slow variables which are related to the streaming 
velocities governing the reversible dynamics. In the present case, the set of slow variables is 
given by {p, g, n^} with the effective Hamiltonian, Eq. ( p^.l] ). The Poisson brackets involving 
the defect density are taken to have the same structure as those involving p, since is 
a scalar quantity. The calculation of the Langevin equations in this case from the effective 
Hamiltonian Eq. is a straightforward generalization of the wavenumber-independent 

case described in Ref. [0,0. The final expressions for the Langevin equations, however, are 
quite comphcated due to the presence of the spatially-varying coefficients in the effective 
Hamiltonian, which will not be presented here. One can put those equations into the field 
theoretic formalism of Martin-Siggia-Rose [Q, which allows one to systematically calculate 



the nonlinear corrrections to the Gaussian response and correlation functions. The detailed 
calulation of the nonlinear corrections relevant to the glass transition problem also follows 
closely the wavenumber-independent case and the reader is referred to Ref. 

In the framework of MOT, nonlinear couplings between density fluctuations in the dy- 
namic viscosity are responsible for the dynamical slowing down. This mechanism is most 
easily seen through the normalized density auto-correlation function: 



0(q, t) = {Sp{q, t)Spi-q, 0))/5(g). (2.6) 

The Laplace transform of 0(q, t), defined by (/)(q, z) = — i dt exp(i2;t)(/)(q, t), is represented 
as 

, / ^ - p^ + ig^r(q, z) 

^^"^^'^ p[z^-Ql{q)]+^zq^T{q^zy ^''^^ 

where ^"^{q) is the microscopic phonon frequency and r(q, 2) is the renormalized dynamic 
viscosity. 

As discussed in Sec. |I|, MCT predicts that <^(q, t), as a function of time, exhibits a 
very elaborate relaxation sequence. As a first step, we restrict the present analysis to the 
intermediate time regime where the MCT analysis is significantly simplified in the following 
aspects. As shown in Ref. |15|, all the structural details of the system in the intermediate 



time regime can be described by a single number A which is related to the power-law and 



the von-Schweidler exponents a and h in the usual way ||T2 

ni-a) _,_ T\l + h) 

r(l-2a) " r(l + 26)' ^ 

where F denotes the gamma function. Consequently, a and h are wavenumber indepen- 
dent, while all the other parameters describing the time relaxation sequence depend on 
the wavenumber. Also, in the intermediate time regime, the defect auto-correlation func- 
tion can be regarded as a constant in time due to its extremely long time scale pf: 
((5?2(q, t)(5n(— q, 0)) = ksT /h"{nct). Thus, the analysis of the system in the a-relaxation 
regime must include the time dependence of the defect auto-correlation function as was ana- 
lyzed in wavenumber independent case in Ref. 0. In the intermediate time regime, this kind 
of complication does not occur and the renormalized viscosity can be expressed in terms of 
the density auto-correlation function alone. 

For the calculation of the nonlinear contribution to F(q, z), we follow the field-theoretic 
calculation in Ref. 0. We find, in the time regime before the primary a— relaxation, the 
dynamic viscosity is represented by 
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r(q,^) 



Tq+ dt e 
Jo 



izt 



_fk_ 

(q, k)0(k, t) + (q, k)0(k, t)0(q - k, t) 



(2.9) 



where V^^^ and V^"^^ are appropriate vertices to be evaluated and To is the bare viscosity 
governing the microscopic dynamics. The appearance of the term hnear in 0(k, t) in Eq. ( p.9|) 
is a direct consequence of the constant defect auto-correlation function in this time regime 



The vertices V^^^ and V^'^^ depend on the coupling constant B and the parameters fi, v 
describing the metastable potential hins) for defects. The basic picture obtained in Ref. |I| 
is that the slowing down corresponds to the weak coupling [B 0) and the low activation 
barrier (yU — > 0) limit with the condition B ~ [j? . In this limit, the vertices are completely 
described by the average number density no and the following two dimensionless parameters: 

no ksT 

where x = {BpQnd)/{nokBT) is the dimensionless parameter describing the scale of the 
coupling energy. We recall that u is the parameter proportional to the third derivative of 
h{nd) at rid = nd- Thus, when the second derivative /x — > 0, the parameter y just represents 
the overall scale of the double well potential. The parameter k is the appropriate ratio 
between the weak coupling (x) and the low barrier (yU^) for defects. We find that for the 
effective Hamiltonian, Eq. the vertices are given by 



s(qj 



y«(q,k) 

no 



no 



where 



f/(q,k) 



2y-U{q, k)s(k)s(q - k) + /€noC(k)s(k) + 
q J 

-^f/2(q, k) - y-U{q, k)l s{k)s{q - k) + 



(q ■ k)noC(k) + (q ■ (q - k))noC(q - k) 



(2.11a) 
(2.11b) 

(2.12) 



We note that in the previous MCT studies of hard-sphere fluids [j^,|T6|-|l^] only the quadratic 
contribution (V^'^^) with the first term in Eq. ( |2.11bD was considered, which corresponds to 
the case where y = k = in this model. 



III. ANALYSIS OF THE MODEL 



Eqs. (|2.7| ) and (2^) with the vertices Eq. ( 2.11 ) completely specify our model which 
incorporates the wavenumber dependence of the system. In this section, we apply the model 
to the case of a hard-sphere fluid using a known approximate but realistic hard-sphere 
structure factor. 



A. Basic Equations from General MCT Analysis 



We first discuss the relevant results from the detailed MCT treatment [|I5|] of Eqs. 
and (|2.9| ) for a general set of vertices V^'^\i = 1,2, ...M. It was found that the microscopic 
structural details of the system are summarized into a single parameter A. This parameter 



is related to the the power-law and von-Schweidler exponents a > and b > via Eq.(2.S). 
The MCT analysis of the coupled set of equations, ( |2.7|) and (|2.9| ) indicates an ergodic- 
nonergodic type transition as parameters describing the vertices approach critical values. 
The nonergodic phase is characterized by 0(q, t oo) = /(q) with /(q) > 0. In the 
present case where M = 2, the nonergodicity parameter /(q) is given as a nonvanishing 
solution of 

7"^= / 7f^[^^'Hq,k)/(k) + V(2)(q,k)/(k)/(q-k)], (3.1) 

i-/(q) {2Try^ J 

while in the ergodic phase, only the trivial solution f{q) = to Eq. (|3.1|) can be found. The 
exponent parameter A is determined for the critical values of the parameters by 

A = /^/ ^e,(q)l^(^)(q,k)(l-/,(k))2(l-/e(q-k))V(k)e,(q-k), (3.2) 

where ec(k) and ec(k) are left and right eigenvectors of 

l^(q,k) = (1 - /,(k))2[v^«(q,k) + 2\/(2)(q,k)/,(q-k)] (3.3) 
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with the normahzation 1 = J d^k/{2TTfec{k)ec{'k) = / ci^k/(27r)3ec(k)e^(k)(l - /^(k)). 

As mentioned in the previous section, MCT predicts a two-step relaxation process in 
the intermediate time regime. After the microscopic time scale tq, the normalized density 
auto-correlation function 0(q, t) shows the power-law decay which can be described as |T7|J2^ 



{ci,t) = fM + hiq)cJ-Y\ (3.5) 



Tn 



This is valid for the time regime tq ^ t ^ Tq, where the characteristic time scale is given by 

Ta ~ |e|~2^. (3.6) 
Here e is the control parameter indicating how far the system is from the transition and 



thus the ideal glass transition corresponds to e = 0. In Eq. (|3.5|) , the coefficients are given 
by Q = -^/jej and h{q) = (1 — fc{.<i)Y(^c{.(l)- For longer times t Ta, the system undergoes 
the von-Schweidler relaxation if e < 

Hci,t) = fM - h{q)A,{-y (3.7) 

with a positive constant This mechanism is valid for the time scale Ta t <^ before 
the a-relaxation, where the characteristic time scale for the a-relaxation is given by 

~ 7 = ^ + 1;. (3.8) 



For e > 0, 0(q, t) decays to a constant, f{q) = fdq) + h{q)c^/y/l — A. We note that the 
critical nonergodicity parameter fdq) and the coefficient h{q) do not depend on the control 
parameter e. 

B. Analysis for a Hard-sphere Fluid 

We calculate the vertices given by Eq. ( |2.11| ), using the Perkus-Yevick solution |^ for 



the hard-sphere structure factor, s{qa) with the Verlet- Weiss correction where a is the 



sphere diameter. For this system, it is more convenient to use the packing fraction 77 instead 
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of the average number density no, where t] = Ticr^no/G. The wavenumber integration is done 
with the mesh size = 300 and the upper cutoff Act = 50. 

As can be seen from Eq. ( |2.11| ), our model for the dynamics of hard-sphere fluids depends 



on the packing fraction rj and the metastability parameters y and k. We recall that y 
describes the overall scale of the double-well potential for defects and k, is the appropriate 
ratio between the weak coupling and the low activation barrier as defined in Eq. (|2.10|) . We 



first investigate whether the double-well potential and the coupling between the defects and 
the density fluctuations can be arranged to yield the slowing down of the dynamics. In order 
to see this, we solve Eq. iteratively for given y, k and rj. We find that for given packing 
fraction 77, a nontrivial solution f{q) > can be found and thus the dynamical slowing 
down is obtained for certain values of parameters y and k. In fact, there exsists a critical 
surface in the y-K space that separates the nonergodic (/(g) > 0) and the ergodic (/(g) = 0) 
phases. This critical surface can be found for various values of the packing fraction (See 
Fig. 1). This result is quite different from that found in the previous MCT studies of the 
hard-sphere system |]^,|TB|-|IB| , which corresponds to the case where y = k = 0. According to 



those studies, the transition occurs at the critical packing fraction rjc ~ 0.525 |TB[. This value 



is determined by the packing fraction at which the first nonvanishing solution of Eq. 
appears when 77 is increased. The present analysis shows that the dynamical slowing down is 
controlled by the metastability parameters y and k, arranging themselves to be on the critical 
surface so that the transition occurs over a range of 77 which contains rjc of the conventional 
MCT. 

For given packing fraction, the values of k on the critical surface do not change very much 
as shown in Fig. 1. As the packing fraction decreases, these critical values of k, increase. Since 
the hard-sphere system remains a fluid for small enough packing fraction, the parameter k 
should be bounded from above. The exact upper bound, however, is not determined in 
the present analysis for reasons discussed below. For high packing fractions rj ^ 0.525, the 
critical surface does not extend to the small y region: For rj ^ 0.525 and for y smaller 
than some value ymm, one always has a nonvanishing solution /(g) to Eq. (|3.1|) . The value 
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of UraiTi increases with increasing 77 (See Fig. 1). We note that for small packing fractions 
{rj ^ 0.505) or for large y {y ^ 1-5), the contribution to /(g) from the large wavenumbers 
become significant as seen in Fig. 2 (a) and (b). This makes the present iteration analysis of 
Eq. ( |3.1| ) with the numerical large wavenumber cutoff Act = 50 less reliable for sufficiently 
small or large packing fractions. However, for the range of packing fractions, 0.51 ^ f] ^ 0.54, 
the present analysis clearly shows the slowing down of the dynamics. 

For the parameters on the critical surface, one can determine the critical nonergodicity 
parameter fdq) and the coefficient h{q) as defined in Eqs. ( p. 51) and (|3.7| ). In Figs. 2 and 



3, /c(g) and h{q) are shown for various values of rj and y along with the corresponding 
quantities of the conventional MCT at 77 = 7]^. In the framework of MCT, approaching 
the critical values of the parameters is represented by the parameter e — > (See Eqs. ( p.5|) 
and ( p.7|) ). The conventional MCT analysis of hard-sphere fluids ||7|,p!^-|T8| is based on the 
assumption that the transition is actually controlled by the packing fraction approaching 
its critical value: e ~ r/ — r/^. Thus, in those analyses, the critical nonergodicity parameter 
fc{(l) and the coefficient h{q) are regarded as parameters independent of rj. In the present 
model, however, the transition occurs when the system becomes metastable with the low 
activation barrier (/i 0) and the weak coupling (x 0) satisfying the condition x/fi'^^K 
(See Eq. (p.lO|) ), while the metastability parameters y and k, arranging themselves to be on 
the critical surface in y-n space shown in Fig. 1. Thus, the control parameter in this case is 
proportional to fi. Therefore, in contrast to the conventional MCT, fdq) and h{q) depend 
on the parameter y describing the overall scale of the defect potential as well as the packing 
fraction. As one can see from Fig. 2 for fixed y, fdQ) shows a variation with 77 such that it 
has larger (smaller) values for low (high) wavenumbers as the packing fraction increases. 

For given point on the critical surface, {y, K{y)), one can calculate the exponent parameter 
A using Eq. ( p.2| ). Thus, in this model A is determined as a function of r] and y as shown 
in Fig. 4. Consequently, the exponents a and b characterizing the intermediate time regime 
are not subject to a particular value of packing fraction, instead they are given as smooth 
functions of 77 with one additional parameter y (See Fig. 5). As a function of the packing 
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fraction, A (77) generally takes smaller values as one goes into the region of either small or large 
packing fractions. It should be noted in Fig. 4 that the values of A become less reliable for 
large y {y ^ 1.5) and small rj {rj 0.505) because of the reason discussed above concerning 
the contribution from the large wavenumbers. It is important to note that within the band 
of packing fractions 0.51 ^ i] ^ 0.54 where the dynamical slowing down is obtained, A(?7) is 
not a constant as a function of rj. 



IV. DISCUSSION AND CONCLUSION 

In this paper, we have been able to show how the metastability parameters describing 
the defects are incorporated in the structural arrest of hard-sphere fluids. The basic picture 
we obtain is that the transition is not a sharp one at a special value of packing fraction, but 
rather the transition occurs over a range of packing fractions as the metastability parameters 
arranging themselves to be on the critical surface as shown in Fig. 1. 

In some sense, the present model just replaces the role played by the critical packing 
fraction rjc ^ 0.525 of the conventional MCT of the hard-sphere system p|JT6|-p^ in the 



slowing down of the dynamics by a band of packing fractions, 0.51 ^ rj ^ 0.54, which include 
?7c. However, our model has very different physical implications from the one studied in the 
conventional MCT in the following aspects: (1) In the conventional MCT, the exponent 
parameter A is defined only at the critical packing fraction rj^, where its value is evaluated 
as A ^ 0.758 Jl^. This value is somewhat larger than most of the values of A found in the 
present model (See Fig. 4). Correspondingly, the values of the power-law and von-Schweidler 
exponents a and b in this model are larger than the conventional MCT values, a = 0.301 and 
b = 0.545 More importantly, since our model is not tied to the notion of rjc, the exponent 
parameter A can be defined over a range of rj. Furthermore, as a function of 1], A is not a 
constant, instead, it shows a variation as shown in Fig. 4. (2) Other physical parameters 
such as the critical nonergodicity parameter fdq) and the coefficient h{q) describing the 
time relaxation sequence also depend on r] in our model, while fdq) and h{q) are defined 
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only aX 7] = Tjc in the conventional MCT. These differences can be summarized in terms of 
the control parameter e of each model such that in the conventional MCT, e ~ r/ — r/^ with 
the well-defined rjc and e is proportional to the metastability parameter /i in the present 
model. 

Recently, there have been many experiments on colloidal systems consisting of spherical 
particles |p6H28||, to which the MCT analysis of hard-sphere fluids has been applied. It is 
claimed in those experiments that the critical packing fraction rjc for the glass transition 
could be identified. The obtained values of rjc range from 0.555 ~ 0.575 for suspensions of 
PMMA particles p6| , |27|l to 0.636 for microgels ||28|. However, it is important to note that 



the identification of rj^. in those experiments is carried out using somewhat indirect methods: 
For the PMMA system, rjc is taken as the value of r] at which a homogeneously nucleated 
crystallization does not occur. For the system studied in Ref. p8[, rjc is determined by fitting 



the experimental data against MCT results Eqs. (|3.5| ) and ( p.7|) assuming that the control 
parameter is e ~ r/ — r/c- Thus, these experimental results should not be considered as a 
conclusive evidence for the existence of ?7c, which is defined by the packing fraction at which 
the density auto-correlation function 0(q, t) shows a sharp ergodic-nonergodic transition, 
i.e. 0(q, t) decays to a finite quantity if t] > t]c and to zero otherwise as t ^ oo. As we shall 
show below, one could instead find an evidence from those experimental data that supports 
the basic picture described in this paper where the glass transition in the hard-sphere system 
is in fact not a sharp one at the well-defined transition density, but it results from the system 
becoming metastable. 

As can be seen from the claimed values of rjc in those experiments, the values of packing 
fraction where the dynamical slowing down is observed are generally higher than the band 
of rj, 0.51 ^ 1] ^ 0.54 obtained here and also than rjc — 0.525 of conventional MCT. This 
discrepancy has usually been neglected in the conventional MCT studies of the hard-sphere 
system |T6|-p!8[]. In fact, in those analyses, only the difference rj — rjc is used to compare with 
the experimental data, where rjc of the conventional MCT is identified with the claimed 
values in the experiments. This discrepancy in t] values between MCT and experiments is 
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usually attributed |jT6|] to the fact that MCT, without including the cutoff effect [0, tends 



to overestimate the tendency to freeze. We believe that this explanation also holds in the 
present case where rjc is replaced by the band of packing fractions. We note that the systems 
considered in those experiments do not have the ideal hard-sphere type potential as their 
inter-particle interactions, which might also contribute to this discrepancy. 

In order to test experimentally the present picture of the dynamical slowing down in 
the hard-sphere system, one must study the dependence of the exponent parameter A and 
the critical nonergodicity parameter fdq) on the packing fraction 77. A variation of these 
quantities as a function of 77 would imply that the system undergoes a smooth transition 
without the special transition packing fraction. In the experiments mentioned above, A 
and fcil) are obtained by using Eqs. ( p.5|) and (|3.7] ) to fit the intermediate time relaxation 



data of 0(q, if:). But, unfortunately, in most cases p6| , p7| , the analysis is carried out under 
the assumption that e ~ 77 — ?7c. In particular, the value of A is not obtained from the 
experimental data, instead the conventional MCT value A ~ 0.758 is directly applied to 
the analysis of the data. Also the apparent variations of fc{(l) with rj observed in Ref. 
P^ , |27[] are not systematically studied but neglected in the analysis. Thus, in order to see 
the 77-dependence of A and fdq), a more careful study of the experimental data is needed 
that treats A, fc{q), h{q)ce and in Eqs. ( ^75|) and ( p.7|) as adjustable parameters without 
assuming any particular form for the control parameter e. 



In fact, this kind of analysis has been performed at one point of the discussion in Ref. [28 



For fixed wavenumeber q, one can indeed see variations of A and fdl) values as functions 



of 7] (See Fig. 5 of Ref. p8[) when they are treated as free parameters. We note that the 
amount of variation of the exponent values is actually comparable to the present findings 
(See Figs. 4 and 5), although there are some differences between the actual values of the 
exponents a and b and the range of packing fractions over which the exponents are measured. 
The authors of Ref. |^ actually concluded that these variations are not systematic ones 
and treated A and /c(q') as constants throughout their analysis. In order to see if this kind 
of data really shows a dependence of the physical parameters, A and fc{q) on t] as described 
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in the present analysis, which results from the metastable nature of the transition, a more 
detailed study is needed covering a wide range of values of the wavenumbers and packing 
fractions. 

The present model predicts that fdq) and A depend on rj with one additional parameter 
y (See Figs. 2 and 4). We note that y is the dimensionless parameter which represents 
the overall scale of the potential energy for the defects. We interpret this metastability 
parameter y as something that depends on the microscopic details of the system and that 
may also have the packing fraction dependence in it, which is not, however, determined 
within the present theory. We believe that the parameter y could be specified in principle as 
a function of packing fraction through a comparison with more detailed experimental data 
showing the dependence of fcig) and A on rj. 

We emphasize that in this paper we are mainly concerned with the qualitative picture 
of the effect of the defect structure on the slowing down of the dynamics of the hard-sphere 
systems, since the defect structure has been approximated in a simplest way using Eq. (p. 51). 
Nevertheless, our model provides an important picture of smooth transition with the physical 
parameters such as the exponents a and b governing the time relaxation sequence and the 
critical nonergodicity parameter fdQ) depending smoothly on the packing fraction. And 
this general picture might be robust in more complex systems. An example of such system 
where the exponents depend smoothly on the temperature can be found in the experiment 
in Ref. §. 

We also note that we have focused on the intermediate time regime before the primary 
a-relaxation. In this time regime, the defect auto-correlation function can be regarded as 
a constant in time |jl|J^ due to its long time scale. For the wavenumber independent case, 
one could extend the FNH formalism by considering the time dependence of the defect 
auto-correlation function as shown in Ref. [Q. One can in principle extend the present 
wavenumber-dependent model to the primary relaxation regime by considering the full time 
and wavenumber dependences of the density and the defect auto-correlation functions fol- 
lowing Ref. [^. The analysis of such model, however, will be a very difficult task. 
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FIGURES 

FIG. 1. The critical surface in the y — k space for various values of packing fraction rj. The 
thick dashed line is ymin(^)- For rj > 0.53, the critical surface does not extend to the small y 
{y < ymin) regions. 

FIG. 2. (a) The critical nonergodicity parameter /c(q) for y = 0.1 and for various values of 
packing fraction (pf). (b) The critical nonergodicity parameter fdQ) for packing fraction r] = 0.52 
and for various values of y. The thick solid line is fcig) of the conventional MCT at rj = rjc — 0.525. 

FIG. 3. The coefficient h{q) describing the power-law and the von-Schweidler relaxations as 
expressed in Eqs. ( |3.5D and ( |3.7D for y = 0.1 and for various values of packing fraction (pf). The 
thick solid line is h{q) of the conventional MCT at r] = r]c ~ 0.525. 

FIG. 4. The exponent parameter A as a function of packing fraction rj for various values of y. 
The star (*) represents the conventional MCT value A ~ 0.758 at rj = tJc — 0.525. 

FIG. 5. The power-law and von-Schweidler exponents a and b as functions of packing fraction 
rj for various values of y. 
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